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Abstract A generating function is given for the number, E(l,k), of irreducible k- 
fold Euler sums, with all possible alternations of sign, and exponents summing to I. Its 
form is remarkably simple: Yl n E(k + 2n,k)x n = J2d\k M^) (1 — x d )~ k ' d /k, where /i is 
the Mobius function. Equivalently, the size of the search space in which /c-fold Euler 
sums of level I are reducible to rational linear combinations of irreducible basis terms is 
S(l, k) = J2 n<k (L(«-H»-i)/2J). Analytical methods, using Tony Hearn's REDUCE, achieve 
this reduction for the 3698 convergent double Euler sums with I < 44; numerical methods, 
using David Bailey's MPPSLQ, achieve it for the 1457 convergent fc-fold sums with / < 7; 
combined methods yield bases for all remaining search spaces with S(l, k) < 34. These 
findings confirm expectations based on Dirk Kreimer's connection of knot theory with 
quantum field theory. The occurrence in perturbative quantum electrodynamics of all 
12 irreducible Euler sums with I < 7 is demonstrated. It is suggested that no further 
transcendental occurs in the four-loop contributions to the electron's magnetic moment. 
Irreducible Euler sums are found to occur in explicit analytical results, for counterterms 
with up to 13 loops, yielding transcendental knot-numbers, up to 23 crossings. 



1 Introduction 



Recent progress in number theory [p], |2| interacts strongly with the connection between 
knot theory and quantum field theory, discovered by Dirk Kreimer p, and intensively 
investigated to 7 loops || by analytical and numerical techniques. The sequence of irre- 
ducible non-alternating double Euler sums studied in [l|] starts with a level-8 sum that 
occurs in the 6- loop renormalization of quantum field theory ||, where its appearance is 
related || to the uniquely positive 8-crossing knot 8i 9 ; the sequence of irreducible non- 
alternating triple Euler sums in || starts with a level- 11 sum that occurs at 7 loops, where 
its appearance is associated with the uniquely positive hyperbolic 11-crossing knot ||. 

These exciting connections, between number theory, knot theory, and quantum field 
theory, led to work with Bob Delbourgo and Dirk Kreimer || , on patterns 0] of transcen- 
dentals in perturbative expansions of field theories with local gauge invariance, and with 
John Gracey and Dirk Kreimer ||, on transcendentals generated by all-order M results 



in field theory, obtainable in the limit of a large number, N, of interacting fields [[U], [11 

In the course of the large- iV analysis [EJ, the number theory in pD] appeared to consti- 
tute a severe obstacle to the development of the connection between knot theory and field 
theory. From the skeining of link diagrams that encode the flow of momenta in Feynman 
diagrams, we repeatedly obtained a family of knots, associated with the occurrence of 
irreducible Euler sums in counterterms. The obstacle was created by the (indubitably 
correct) 'rule of 3' discovered in [|l[], for non-alternating sums ((a,b) = X] n>m l/fi'm* of 
level I = a + b. The analysis of |l] shows that non-alternating sums of odd levels are re- 
ducible, while at even level / = 2p+2 there are [p/3j irreducibles, where |_- • -J is the integer 
part. At levels 8 and 10 this made us very happy, since we had the 8-crossing knot 819 
to associate with £(5,3), and the 10-crossing knot IO124 to associate with £(7,3). There- 
after the knots increase in number by a 'rule of two', giving \_p/2\ knots with / = 2p + 4 
crossings. So there are two 12-crossing knots, while [I| has only one level-12 irreducible, 
and two 14-crossing knots, which is [|l| the number of level- 14 irreducibles. 

Faced with a 12-crossing knot in search of a number, we saw two ways to turn: to 
study 4-fold non-alternating sums, or 2-fold sums with alternating signs. The first route 
is numerically intensive: it soon emerges that well over 100 significant figures are needed 
to find integer relations between 4- fold sums at level 12. The second route is analytically 
challenging; it soon emerges that at all even levels I > 6 there are relations between 
alternating double sums that cannot be derived from any of the identities given in |L|. 

Remarkably, these two routes lead, eventually, to the same answer. The extra 12- 
crossing knot is indeed associated with the existence of a 4-fold non-alternating sum, 
£(4,4,2,2) = £ 

n>m>p>q^/ n ' i ' mi p 2( l 2 ^ which cannot be reduced to non-alternating sums 
of lower levels. It is, equivalently, associated with the existence of an irreducible alter- 
nating double sum, U 9j3 = En>m {(-l) n /n 9 } {(-l) m /m 3 }. The equivalence stems from 
the unsuspected circumstance that the combination £(4,4,2,2) — (8/3) 3 f/ 9i3 , and only 
this combination, is reducible to non-alternating double sums. Moreover, Z = 12 is the 
lowest level at which the reduction of non-alternating 4-fold sums necessarily entails an 
alternating double sum. The 'problem pair' of knots are a problem no more. Their entries 
in the knot-to-number dictionary |12] record that they led to a new discovery in number 
theory: the reduction of non- alternating sums necessarily entails alternating sums. 
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This discovery led me to study the whole universe of fc-fold Euler sums, with all 
possible alternations of sign, at all levels I. As will be shown in this paper, it is governed 
by beautifully simple rules, which might have remained hidden, were it not for Dirk 
Kreimer's persistent transformation of the Feynman diagrams of field theory to produce 
a pair of 12-crossing knots. 

The remainder of the paper is organized as follows. Section 2 statesQ the formula for 
the number, E(l,k), of irreducible fc-fold sums, with all possible alternations of sign, at 
level I. Section 3 outlines the process by which it was discovered. The anterior numerics 
of Section 4 describe the high-precision evaluation methods and integer-relation searches 
that helped to produce the formula; the posterior analytics of Section 5 describe computer- 
algebra proofs of rigorous upper bounds onE(l, k) that are respected (and often saturated) 
by it. Section 6 summarizes numerical and analytical findings by listing an instructive 
choice of concrete bases. Section 7 shows that all 12 of the irreducible sums with I < 7 
appear in perturbative quantum electrodynamics. Section 8 considers the import, for 
quantum field theory, for knot theory, and for number theory, of results obtained by 
calculations up to level 23, corresponding to knots with up to 23 crossings, and to Feynman 
diagrams with up to 13 loops. 



2 Result 

To specify an alternating fc-fold Euler sum, one may give a string of k signs and a string 
of k positive integers. It is veryQ convenient to combine these strings, by defining 

CK . . . , a k ) = £ 11 ^ — , (1) 

ni>n i+ i i=l Tl i 

for k > 1, on the strict understanding that the arguments are non-zero integers, and that 
a\ ^ 1, to prevent a divergence. Hence one avoids a proliferation of disparate symbols for 
the 2 k types of fc-fold sum. The correspondence with the double-sum notations of [0] is 

CM) = Mt,s), (2) 

C(-s,t) = a h (t,s), (3) 

C{s,-t) = -<r a (M), (4) 

U Stt = C(-8,-t) = -a a (t,s), (5) 

for positive integers s and t, with emphatically no implication of analytic continuation. 



In perturbative quantum field theory, three-loop radiative corrections fl3| , |i~4| involve 

I i \n+m 

U 3 ,= £ LJ|L_ = | C (4)-2{Li 4 (|) + iln 2 2(ln 2 2- 7 r 2 )}. (6) 



n> m >0 n7Jl 



At six [|TJ, [Ti| [Tjj] and seven |J loops, counterterms associated with the (4, 3) and (5, 3) 
torus knots, 819 Q and IO124 0, involve and Ur,3, whose irreducibility is equiva- 
lent to that of C(5,3) and £(7,3), respectively. In higher counterterms, the independent 
irreducibility of Uq^ and Ur,5 is associated with a pair of 12-crossing knots. 

1 It would be very difficult to prove. One cannot even prove that C 2 (3)/C(6) is irrational. 
2 Without such a convention, results such as (E2lB0f) become almost unreadable. 
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In what follows, the number of summations, k, is referred to as the depth of the 
sum (HD, and I = J2i I a i\ is referred to as its level. The set Si^ of convergent sums of level 
I and depth k has 



N(l,k) 



I - 1 



fc -iM-2 



\k-lj \k-2J 
elements, as is easily proven by induction. The number of convergent sums at level I is 

i 

N(l) = k) = 4 x 3 l ~ 2 , ({ 

k=l 

provided that I > 1. 

Table 1: Euler's triangleR of irreducibles, at level / and depth k, for I + k < 32 . 
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3 
7 

14 

25 

42 

66 

99 
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273 



1 

3 
9 

20 

42 

75 

132 

212 

333 

497 



1 

4 

12 

30 

66 

132 

245 

429 

715 



1 

4 

15 
40 
99 
212 
429 
800 



1 

5 

18 

55 

143 

333 

715 



1 

5 

22 
70 
200 
497 



1 

6 

26 
91 
273 



1 

6 

30 
112 



1 

7 

35 



1 
7 



3 It seems appropriate to call it Euler's triangle; see Section 8.3 for the connection with Pascal's. 
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The main purpose of this work is to determine, as in Table 1, the number, E(l, k), of 
irreducible sums in S^k, i-e. the minimum number of sums which, together with sums of 
lesser depth and products of sums of lower level, furnish a basis for expressing the elements 
of S^fc as linear combinations of terms, with rational coefficients. The conclusion is that 

2 /l±k\ 
E(l, k) = 6 l>2 5 k>1 + 7-t E Kd) £ , (9) 

1 + K 2d\l±k \ 2d / 

where the summation is over the positive integers d such that (I ±k)/2d are integers, and 
is weighted by the Mobius function, //(d), which vanishes if d is divisible by the square of 
a prime and otherwise is ±1, according as whether d has an even or odd number of prime 
divisors. With the exception of In 2, from Si i, and 7r 2 , from S^i, the irreducibles come 
from Sk+2j,k, with j > 0. Moreover, Sk+2j,k contains the same number of irreducibles as 
Sj+2k,j, as illustrated in Table 1. 

From (H) one obtains the number of irreducibles at level /: 

E{1) = E E(l, k) = W Kl/d) {F d+1 + Fa-x} , (10) 

k=l 1 d\l 



in terms of the Fibonacci numbers. The integer sequence (pT0|) is tabulated as M0317 by 
Sloane and Plouffe flTSfl , who record its origin in the study of congruence identities. 
It shows that the irreducibles are an exponentially decreasing fraction of the number of 
convergent sums at level I: 

This relative sparsity is illustrated in Table 2. 

Table 2: The numbers, E{1) and N(l), of irreducibles and sums, at level /, for I < 13 . 
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E(l) 
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2 2 


4 


5 


8 11 


18 


25 


40 


N(l) 


1 


4 


12 


36 


108 324 


972 


2916 


8748 26244 


78732 


236196 


708588 



From (0) one may obtain the size, S(l, k), of the search space for sums of level / and 
depth k, i.e. the minimum number of terms that allow one to express every element of Si^ 
as a linear combination, with rational coefficients. This basis consists of the irreducibles 
with level / and depths no greater than k, together with all the independent terms that 
are formed from products of sums whose levels sum to / and whose depths sum to no 
more than k. 

One may generate S(l, k) from E(l, k), by expanding 



in 



E 

?1=1 



E [{X^x^y+j^Y, X lAi x l y h 

1=1 \ k=l i=l 



12) 



to order / max in both x and y, where Z max is greatest level required, and Xi^,i serves as 
a symbol for the ith irreducible in Si t k, so that (-X^i)" stands for any rational multiple 



5 



of 7r 2n , and hence for ((2n). Then S(l, k) is obtained by selecting the terms of order x l , 
dropping powers of y higher than y k , setting y = 1, and counting the length of the resulting 
expression. This procedure is easily implemented in REDUCE, Maple, Mathematica, etc. 

Table 3: The size, S(l, k), of the search space for sums in S^, for k < I < 15 . 
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15 


25 


40 


46 


53 


54 


55 














10 


1 


6 


16 


36 


51 


72 


79 


87 


88 


89 












11 


1 


6 


21 


41 


76 


97 


125 


133 


142 
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144 










12 


1 


7 


22 


57 


92 


148 


176 


212 


221 


231 


232 


233 








13 


1 


7 


28 


63 


133 


189 


273 


309 


354 


364 


375 


376 


377 
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29 


85 


155 


281 


365 


485 


530 


585 


596 


608 


609 


610 




15 


1 


8 


36 


92 


218 


344 


554 


674 


839 


894 


960 


972 


985 


986 


987 



Inspection of Table 3 reveals that the sizes satisfy the recurrence relation 
S(l + 1, k) = S(l, k - 1) + S(l - 1, k) , for I > k > 1 . 
Moreover, the Fibonacci numbers appear at maximum depth: 

S(l,l) = S(l,l-l) + l = F l+1 . 
From the Lucas relation between Fibonacci and binomial numbers, one obtains 



fc-i / 

s(i,k) = J2[ 

n=0 V 




(13) 



(14) 



(15) 



as the solution to ([T^,|T4]). 

Since the process of generating Table 3 from Table 1 is reversible, formula (|9|) may be 
replaced by the simple statement that at level I the size of the search space increases by 
the binomial coefficient 



S(l,k + 1) - S(l,k) 



i+k-i 

2 
k 



(16) 



when the depth is increased from k to k + 1. This pragmatic formulation^ is rather helpful, 
when using an integer-relation search algorithm, such as PSLQ p0| . An even simpler, 
though informal, restatement appears in Section 8.3. 

4 Computation leaves no doubt as to the equivalence of M) and (|l6|), though it is not yet proven. 
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3 Discovery 



Some brief historical remarks seem in order at this stage, since my route to @ in fact 
began with the observation of 7 members of the Fibonacci sequence (|l4l), during the 
course of 800-significant-figure integer-relation searches entailed by the relation between 
knot theory and field theory. 

Using David Bailey's magnificent MPPSLQ implementation of PSLQ p0| , I suc- 
ceeded in reducing all 1457 sums in {S 1 ^ | 1 < k < I < 7} to just 12 numbers and their 
products. These irreducible numbers may conveniently be taken as 

In 2, vr 2 , { C(0 | I = 3, 5, 7}, { a(l) | I = 4, 5, 6, 7}, U 5A , C(5, 1, -1), C(3, 3,-1), (17) 

with the polylogarithms 



l\ { 12 Vln2/ J ^2 n n l v ' 

populating the deepest diagonal of Table 1. As in (|), the definition ( |i~8"D postpones the 
appearance of 7r 2 (ln2)' -2 and (In I) 1 to depths / — 1 and I, respectively, which is required 
by (|i~6|). Moreover ct(l) = a(2) = 0, and a(3) = |C(3), which is not a new irreducible. 

After noting the Fibonacci sequence for the maximum sizes S(l, I), with I < 7, I sought 
a combinatoric form for S(l,k). Formula flT5p suggested itself on the empirical basis of 
the 28 cases with 1 < k < I < 7, and was then submitted to intense numerical and 
analytical tests at higher levels, as indicated in Sections 4 and 5. I implemented the 
generator flT2|), using the weight and length commands of REDUCE |52fl , and constructed 
Table 1, working backwards from ([H|). Next came the observation that the generating 
functions for the k = 2 and k = 3 columns of Table 1 have comparable forms: 

G 2 (x) = \ {1/(1 - xf - 1/(1 - x 2 )} , G 8 (aO = | {1/(1 - xf - 1/(1 - x 3 )} . (19) 

The symmetry of Table 1 was vital to the discovery of the simple formula 

OO 1 

G k (x) = Y,E{k + 2n, k)x n = -Y, Kd) (1 - x d y k / d , (20) 

n=0 fc d\k 

which produces formula (|). Computation of J2k ^{1, k), for 3 < / < 100, revealed that it 
produced the integers nearest to J2d\i ^(l/d)(fi d /l, where = |(\/5+l) is the golden section. 
The equivalent form (|T0|) was obtained by submitting the integer sequence E(3), . . . , E(23) 



to Neil Sloane's helpful on-line version of |18[ . It was noteworthy that this lookup 
returned the values E(l) = E(2) = 1, from the Fibonacci form (|i~0"D , agreeing with the 
appearance of In 2, at I = 1, and it 2 , at / = 2. Their inclusion in the k = 1 column of 
Table 1, above Euler's triangle, thus became appropriate, as they appear to usher in the 
higher transcendentals, in much the same way that F\ = Fi = 1 seed the exponential 
growth of MPPSLQ's CPUtime (or Fibonacci's rabbits) along the deepest diagonal of 
Table 3. 

It remains to describe yet more probing numerical and analytical tests, in further 
support of the claim in (0), and its equivalent version in (|20|). Nonetheless, it is hoped 



that the reader already shares some of my feeling that these two formulae are simply too 
beautiful to be wrong. 
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4 Anterior numerics 



4.1 Numerical evaluation 

Suppose that one wishes to obtain a high-precision approximation to Soo = hnijv^oo S(N), 
where the truncated sum S(N) = J2 n <N R(n) has a summand R{n) with an asymptotic 
series in 1/n, starting at 1/n +1 , with C > 0. From {S(n) \ N — M < n < N}, one may 
form a table {T(n, m) \ N — M + m<n<N, < m < M}, by the procedure 

(C + m + n) T(n, m) - nT(n - 1, m) 
T(n,m + 1) = ^-^ , (21) 

with T(n, 0) = S(n). The method exploits the vanishing of C + m + n — n/(l — 1/n) +T ™ 



as n — > oo. It takes M(M + l)/2 applications of (|2l|) to produce T(N,M). Provided 
that N » M > 1, and that rounding errors have been controlled, one obtains the 
approximation = T(N,M) + 0(M!/iV M+c '), where the factorial becomes significant 
for iV » M » C > 1. 

This elementary and economical method of accelerated convergence is applicable to 
every Euler sum of the form (|l|) that has no argument equal to unity, in which case 



repeated appeal to the Euler-Maclaurin formula [|24| underwrites the absence of logarithms 
in the expansion of the truncation error, and the outermost summation is of the form 
J2 n {R(2n) ± R(2n + 1)}, with R{n) = 0(l/n a ), for an argument ai = ±a. Thus one 
should store truncations after even increments of n\, and set C = a\ — 1 for a\ > 2, or 
C = | ai| for ai < —1, in procedure (0). To obtain the starting values, one has merely 
to set up a single loop that updates and stores each layer of the nest as its particular 
summation variable, rij, assumes the even and odd values within the loop. Thus the 
evaluation time for a truncation at iV of a fc-fold sum is roughly proportional to kN. 

For sums with no unit arguments, one needs therefore only a few lines of conventional 
FORTRAN, which may be handed over to David Bailey's TRANSMP H utility, to pro- 
duce code that calls his MPFUN [26] multiple-precision subroutines. As a rule of thumb, 



the working accuracy should be somewhat better than the square of the desired output 
accuracy, when using (pi]). When, and only when, rounding errors are so controlled, an 
output accuracy of very roughly M\/N M is achieved by M iterations of (|2T|) , with input 
data obtained from looping over iV pairs of successive even and odd integers. For a 4-fold 
sum, the accumulation of data takes roughly 50iV calls of MPFUN subroutines, and the 
acceleration of the convergence takes roughly 8M 2 calls. Thus, to achieve P significant 
figures for a 4-fold sum, one should choose a value of M that keeps the time factor, 
T w 50(10 P M!) 1 / M + 8M 2 , close to its absolute minimum. I find that between 780 and 
800 significant figures are reliably and efficiently achieved with M = 440 iterations, and 
truncation at N = 10 4 , which entails T«2x 10 6 calls to MPFUN subroutines, operating 
at multiple precision 1700. This takes less than half an hour on a DEC Alpha 3000-600 
machine, corresponding to a call rate that is faster than 1 kHz. The memory requirement 
is less than 1 MB: an array of 440 x 240 4- byte cells holds the truncations in multiple 
precision and is updated iteratively by (|21"|). 

Of course, the above method fails as soon as one sets any of the arguments to unity. 
However, all is not lost. By iteratively applying the Euler-Maclaurin formula, one arrives 
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at the conclusion that the maximum power of Inn in the truncation error YUj Aj(hm)V nJ 
is the largest number, d, of successive units in the string of arguments, since only for 
j = 1 does the integration of (InxY/x^ increase the power of lnx. Consider, for example, 
((a, 1,6, 1), with a / 1 and 6^1. By the time the logarithm generated by the 
summation is felt by the n 2 summation, it has acquired an inverse power of n 2 , thanks to 
the benignity of the Euler-Maclaurin formula for the n 3 summation. Thus this sum has 
(iemon-number d — 1, even though it contains two unit arguments. On the other hand, 
C(a, 6, 1, 1) and ((a, 1, 1, 6) have d — 2, as does ((a, 1, 1, 6, 1, a, 1, 1), for example. At level 
/ the most demonic convergent sum has unity for all its arguments, except for the first, 
which must therefore be a± = — 1. Fortunately, it is possible to give exact expressions for 
this beast, with d — I — 1, and another, with d = 1 — 2. With a string of n unit arguments 
denoted by {l} n , the all-level results 

c(-i,{ih-i) = L 7r ^, -c(-i,-i,{iw = u(±) = £ — , (22) 

1 ■ n=l z n 

were inferred numerically, and then obtained analytically, by Jon Borwein and David 
Bradley, as the coefficients of t l in G(t, —1) and t dzG(t, z)/(l + z), generated by the 
trivially summable hypergeometric series G(t, z) = tz 2 -Fi(l + t, 1; 2; z) = (1 — z)~ l — 1. 
From t Jq dzG(t, z)/z one generates the corollary, £(2, {l}/_ 2 ) = C(0? °f a theorem 

£ C(ai, a 2 , . . . , a k ) = C(0 , (23) 

from Andrew Granville, which was used to check evaluations of non-alternating sums. 

To mitigate the computational difficulties caused by unit arguments, truncation errors 
may be obtained analytically for combinations of Euler sums of the form 

^ (sign a)" ^ (sign 6,)^ 
S(a; 61, ... , 6 fc _i) = £ J] L ■ ( 24 ) 

n=l ' 4 i=l mi=l 

It seems appropriate to call p^ ) a 6oxed sum, since the symmetrical inner summations 
span a lattice that is confined to a (k — l)-dimensional hypercube by the outer summation 
variable. It is built out of symmetrical combinations of Euler sums with depths no greater 
than k, and is, so to speak, a 'cheap boxed set', available at a reduced computational 
price, since it requires only the multiplication of k — 1 polygamma Euler-Maclaurin series, 
followed by a single further application of the Euler-Maclaurin formula, to determine the 
truncation error to any order for which ones favorite computer-algebra engine has the 
power to multiply, differentiate, and integrate double series in Inn and 1/n. Only depth-1 
data is needed, since the inner summations can be rewritten as Ysm- ~ J2 mi >m with the 
first term giving a constant and the second an asymptotic series in 1/n, except for 6j = 1, 
which gives Euler's constant, a log, and an asymptotic series. 

This Euler-Maclaurin method is the obvious generalization of that used in ||28||, in the 



restricted cases with 6j = 1, or 6j = —1. In the general case, much computer-algebra time 
may be consumed by multiplying asymptotic series for the distinct values of 6j and then 
integrating a long expression involving many powers of In n and 1/n, though the algorithm 
is straightforward to implement. To obtain a few hundred significant figures in a short 
time one may stay within REDUCE, without feeding thousands of lines of FORTRAN 
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statements to TRANSMP; at higher precision one benefits from this transfer of function, 
at the expense of non-trivial file management, when handling many boxed sums. 

From an analytical point of view, it is highly significant that every 3-fold Euler sum 
can be transformed into a boxed sum: C( a i> a 2 , 03) — —S(a 2 ; ai, a 3 ), where ~ will hence- 
forward stand for 'equality modulo terms of lesser depth, their products, and consequent 
questions of convergence'. The argument is simple f2|: to the summation with N > n± > 
n 2 > n%> 1, for the truncation of C( a i ; a 2, 03) , add that with N > n 2 > 711,3 > 1, for the 
truncation of S(a 2 ; a±, a 3 ). Adding the case n\ = n 2 , which has lesser depth, one has cov- 
ered the values N > n\ > 1, and hence has a product of truncated sums of lesser depth. 
Thus C(ai,a 2 ,a 3 ) + 5 , (a 2 ;ai,a 3 ) ~ 0. The questions of convergence obviously concern 
the case with a 2 = 1, in the limit N — > 00. Such issues are handled with great dexterity 
in 0, in the case of non-alternating triple sums. Dealing, as now, with truncated sums, no 
problem of convergence arises. The snag is that one must multiply truncations to obtain 
the product term, so the process becomes both messy and ad hoc, from the perspective of 
a programmer seeking a systematic algorithm for sums of any depth. From a numerical 
point of view, little is gained from boxability at depth k = 3. 

Truly nested Euler sums ([I]) have depth k > 4, where only combinations of them 
can be boxed. In the next subsection it will be shown that at k = 4 there first occurs 
a significant phenomenon, of central relevance to the claim of Section 2, and to knot- 
theoretical studies H, H H ||. But first, there is an outstanding computational dilemma 
to confront. Should one attempt to automate the process of chaining applications of 
the Euler-Maclaurin formula, feeding in the required values for multiple Euler sums from 
runs of lesser depth, and carefully separating the odd and even summations, at each 
link of a chain whose length one would like to vary? Or should one, rather, use the 
empirical truncation data to accelerate the convergence? In the absence of logs, the one- 
line procedure ( plj) settles the issue, to my mind. Why use masses of computer algebra 
to feed into TRANSMP information that is already sitting in the numerical data, in such 
immediately usable form? When there are logs, from unit arguments, the matter is moot. 
Having taken the empirical approach when the going was easy, I opted to stick with it, 
when the going got tough, at depths k > 4, with al > demons. 

It is clear that high-precision knowledge of K(d, M) = 1 + (d + 1)M truncated values, 
for a sum with demon-number d, should suffice to accelerate the convergence by a factor 
of roughly M\/N M , as was achieved with (p!~D, for d = 0. How best to achieve this, for 
d > 0, is another matter. Unable to devise an easily programmable iterative method 
like I returned to the brute- force method of using K(d,M) truncations to solve 

directly for K(d, M) unknowns, as was done a decade ago [pl| in the investigations which 
first suggested the appearance of irreducibles with depth k > 1 at the six-loop level of 
renormalization of quantum field theory, a prediction amply confirmed by recent six- and 
seven- loop analysis f| , and illuminated by knot theory || f|] . 



After translating the Gauss- Jordan |29| method into MPFUN calls, one may system- 



atically obtain the K(d } M) coefficients, in the approximation 

d C+M-l a 

S(n) « + J2(lnnY £ ~J~ > ( 2 5) 
,=o J= c nJ 

from K(d,M) truncations, and in particular find an accurate value for S^. I chose the 
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truncations {S(N + pL) | < p < (d + 1)M}, taken distances L apart, starting at some 
value N. The choices of M, N, L, and working accuracy, to achieve a desired output 
accuracy for a sum with demon-number d, are an art learnt by experience, not yet a 
science that is fit to be explained here. Suffice it to say that it is wise to work at the cube 
of the desired accuracy, and that solving 10 3 sets of 10 3 equations, at a working accuracy 
of one part in 10 103 , in order to perform 10 3 Euler-sum searches with the MPPSLQ [ fZT| ] 
integer-relation finder, is demanding of core memory, CPUtime, vigilance, and patience. 
The rewards, in increased analytical understanding, are considerable. 



4.2 Exact numerical results 

By an exact numerical result, I mean an equation whose exactness is beyond intelligent 
doubt, yet is validated, to date, only by very high precision numerical evaluation. It is, 
strictly speaking, possible that an equation presented here is a divine hoax, just as the 
rationality of ( 2 (3)/((6) is still a possibility. The improbability beggars all description. 

The exact numerical result that sparked the genesis of Tables 1 and 3, and revealed 
their utter simplicity, was found at level I = 12 and depth k = 4. It reads 

2 5 • 3 3 C(4, 4, 2, 2) = 2 5 • 3 2 ( 4 (3) + 2 6 • 3 3 • 5 • 13 C(9) C(3) + 2 6 • 3 3 • 7 • 13 C(7) C(5) 
+ 2 7 • 3 5 C(7) C(3) C(2) + 2 6 • 3 5 C 2 (5) C(2) - 2 6 • 3 3 • 5 • 7 C(5) C(4) C(3) 

-2«.3 2C (6)C 2 (3)-^=C(12) 

+ 2 4 • 3 3 • 5 • 7C(6, 2) C(4) - 2 7 • 3 3 C(8, 2) ((2) - 2 6 • 3 2 • ll 2 ((10, 2) 

+ 2 14 f/ 9 , 3 (26) 

and the sting is in its tail. 



Exceptionally, factorizations of rationals were written in (|26|) , with a central dot (not 
to be confused with a decimal point) denoting multiplication. If one dislikes the number 
691, one may remove it, using ((12) = 69l7r 12 /(3 6 -5 3 -7 2 -ll-13). Factorizations were given 
above lest the reader found the unfactorized rationals implausible, on first encountering 
an exact numerical result obtained by MPPSLQ. The practice will not be continued. 

It is apparent that one needs at least 100 significant figures to discover (^6j) , because 
of the product of two 5-digit primes in the 7r 12 term. MPPSLQ (almost always) finds the 
integer relation, Y,i n iSi ~ 0, with smallest euclidean norm, (J^irf) 1 ^ 2 , consistent with 
the requested accuracy of fit. If one knew ((4,4,2,2) to only 100 significant figures, the 
routine would be at perfect liberty to return 13 'random' 8-digit integers that happened 
to fit at 100 significant figures. It would not care that the true form has attractive 
factorizations of all integers save one. At 800 significant figures, which is the accuracy to 
which (EBf) has been validated, the probability of it being spurious is of order 10 -700 . 



The import of ( |2"6] ) is dramatic: non-alternating sums, with exclusively positive argu- 
ments in (0), do not inhabit a cosy little world of their own, uncontaminated by contact 
with their alternating cousins, as the presence of C/g^ clearly demonstrates. 

As explained in the introduction, this was wonderful news for the connection between 
knot theory and quantum field theory 0, [|, || ||. It was also what sparked the present 
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systematic inquiry into the whole universe of Euler sums of the form (JI|), setting firmly 
aside the notion that there is something special about non-alternating sums. The liber- 



ating effect is apparent in the discovery of the simple formulae (|16],^0D, which give order 
to the larger universe so embraced. 

Inspection of Table 3 reveals some practical limits to exploration. Even with 800 
significant figures one might expect to encounter problems finding relations in £15,3, or 
610,4, each of which has 36 basis terms. If just one of the 37 integers in a desired integer 
relation exceeds 10 23 , then MPPSLQ may (quite properly) fit 800-significant-figure data 
with 37 'random' 22-digit integers. Experience shows that integers of order 10 22 are 
produced by successful searches in 613,3, with 28 basis terms. So investigation of 615,3 
was judged to be imprudent, with 'only' 800 significant figures at hand. Because of the 
importance of level 10 in 0, investigation of 6104 is reported in Section 6. 



The successful fit at all levels up to and including I = 7, using the 12 numbers (|Tj 
has been reported, as have the all- level results (^). There remains the apparently trifling 
matter of an exact numerical relation in Se,2 : 

t/4, 2 = IC(6)-|C 2 (3), (27) 

with coefficients that Euler could, no doubt, have found by mental arithmetic. Surpris- 
ingly, nothing in the most recent work on double sums |]] suggested the existence of such 
a relation. It turned out to be relatively easy to devise an analytical proof, when thus 
apprised of its need, so the subject is postponed to the next section, where the armory of 
analytical tools is augmented. Yet the relation belongs here, since it was David Bailey's 
engine that disclosed it. The analytical techniques, developed to derive (p7|), feed back 
useful information for further numerical analysis. 

First, they clear up the whole of the double-sum sector, for good and all, confirming 
the knot-theoretic expectation of a 'rule of 2', with E(2p,2) = [p/2\, and thus allowing 
numerical exploration to progress to A; > 3. Secondly, they provide rigorous (though 
non-optimal) upper bounds: 

£7(7,3) < 3, £7(8,4) < 4, £7(9,3) < 6, £7(11,3) < 11, £7(13,3) < 17. (28) 

Thirdly, for each bound it is possible to find an overcomplete basis, whose size is deter- 
mined by the bound. Thus one needs to evaluate only a small fraction of the sums in 
these sectors, and then use MPPSLQ to reduce an overcomplete basis to a minimal basis 
of size S(l,k). If the formula for S(l,k) were false, MPPSLQ might sometimes fail to 
reduce the basis down to the claimed size, or might reduce it down to a size smaller than 
that claimed by (|T5|). Of course, it never did the latter, else the claim would not have 
been made. The fact that it never did the former, with 800-significant-figure data, in 
search spaces of sizes up to 36, is a testament to its author as well as to the formula. 
Finally, and most fortunately, it is possible to construct overcomplete bases, in the spaces 
bounded by (p8|), that are demon-free. Hence 800 significant figures are available, in less 
than half an hour per sum, at the touch of button (|2l|). 

Thanks to these circumstances, all Euler sums in {6^3 \ I < 14} and {6^4 | I < 9} 
have been shown to be reducible to bases of the sizes given in Table 3, by the operation of 
MPPSLQ on overcomplete bases. Relatively simple examples of such reductions, in 69 3 , 
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611,3 and 5*8,4, are provided by: 

C(3,-3,-3) = 6C(5,-l,-3) + 6C(3,-l,-5)-^ln2C(3)C(5) + 6f/ 5 ,iC(3) 

+ ^ C(2) C(7) - | C 3 (3) + f| C(3) C(6) + If C(4) C(5) - 2 ~W C(9), 

C(3,-5,-3) = i§^C(5,3,3) + 15C(7,-l,-3) + 15C(3,-l,-7) + ^f/ 5 , 3 C(3) 

+ 15 f/ 7il C(3) - If In 2 C(3) C(7) - W C(H) + ^ C(2) C(9) 
" C 2 (3) C(5) + Hi 1 C(3) C(8) + » C(4) C(7) + W C(5) C(6), 

C(3, -1,3,-1) = | C(-3, -3, -1, -1) - f C(-5, -1, -1, -1) - # 1/5,1 C(2) 



fgf [7 5 , 3 + ^ U 7>1 + f ln^2 C/ 5>1 + § ln^ 2 C a (3) 
5 f In 2 2 C(6) - In 2 C(2)C(5) - fUn 2 C(3) C(4) 
W C(2) C 2 (3) + ™p CO) C(5) + fPf C(8). ( 20 ) 



In 5i3,3, the relations are more complex, with the prime factor 102149068537421 ap- 
pearing in one case. Nonetheless, the probability of a spurious fit is less than 10 -200 , 
in all cases, and is often much less than this. The existence of further relations, forbid- 
den by (0), cannot be disproved by numerical methods. The euclidean norms of such 
unwanted relations would, however, exceed those of the discovered relations by factors 
ranging between 10 10 and 10 20 , which makes it extremely implausible that the formula is 
in error in any of the spaces with 5(7, k) < 5(8, 8) = F 9 = 34. 



5 Posterior analytics 
5.1 Analytical tools 

The analysis of |T], |2| makes use of two very simple types of relation between Euler sums. 
In the general case of fc-fold sums, with all possible alternations of sign, it is somewhat 
difficult to notate these relations, in all generality. To avoid cumbersome formulae, terms 
that involve sums of lower depth, and their products, will be omitted, as in the case of 
C(ai,a 2 , a 3 ) ~ -5(a 2 ;ai,a 3 ). 

The first type of relation involves permutations of arguments: 

~ C(«i> 02, o 3 , • • • , a k ) + C(a 2 , ai, a 3 , . . . , a k ) + £(a 2 , a 3 , a 1} . . . , a k ) + . . . 

+ C(°2, a 3 , a 4 , . . . , a x ) . (30) 

The proof is trivial: by including all insertions of a% in the string a 2 , a 3 , . . . , a*,, one obtains 
a combination of sums that differs from the product C(°i) C(°2, 03, ■ ■ ■ , Ofc) only by terms 
in which summation variables are equal, corresponding to sums of depth k — 1. For 4- fold 
sums, such relations reduce a set of 24 possible permutations to a set of 9, when the 
arguments are distinct. 

The second type of relation follows from use of the partial-fraction identity [I], [2|] 
1 v 1 ( fa + b - s - 1\ 1 



A a B h (A + B) a+b - s \ \ a - s J A s 
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for positive integer a and b. To see how this is used, consider the product ((a)((b,c,d), 
with positive arguments. It may be written as J2 n , P , q ,r ^/n a (p + q + r) b (q + r) c r d , where 
each summation variable runs over all the positive integers. Now apply (0), setting 
A = n and B = p + q + r. The second type of resulting partial fraction is of the form 
Y^npqr 1 / ( n + V + Q + r) a+fe_s (p + g + r) s (q + r) c r d , which is an Euler sum. To the first 
type, apply ( pHf ) with A = n and B = q + r. Its second terms are also Euler sums. To 
its first, apply (0) with A = n and B = r. Each term so produced is an Euler sum. 
Thus one has obtained a relation for non-alternating sums. By including signs, 16 such 
relations can be generated. In general, one gets 2 k relations by k — 1 applications of (|3l"l) 
for every set of k exponents that one chooses for the initial product of sums. 

It can seen that there is no scarcity of trivially derivable relations between Euler 
sums. The notable achievement of [0J was to organize the relations between double sums 
in such a way as to prove the reducibility of all double sums of odd level. In non- 
alternating triple sums of even level were proven to be reducible. It was conjectured that 
non-alternating sums of level / and depth k are reducible whenever I + k is odd. The 
stronger claim made by ((^) is that this applies to alternating sums as well. In the course 
of the present work, reducibility has been demonstrated, by a combination of analytical 
and numerical methods, for all odd values of I + k such that S(l, k) < 5(14, 3) = 29. 

As remarked previously, the identities of pfl are insufficient to derive the simple re- 



lation (27). No tally was given in |l| of the numbers of alternating double sums left 
unreduced at even levels, though the tally \j>/3\ was made for non- alternating double 
sums at level 2p + 2. Using REDUCE, one easily discovers that the relations given in 
allow reduction of double sums to the set { U n+ 2 m ,n I min(n, m) > 0} and that no further 
reduction is possible without additional input. Using MPPSLQ, one easily discovers that 
a truly irreducible set is furnished by { U2 a +3,2b+i | a > 6 > 0}. Thus the relations derived 
in [0 are insufficient in a way that is very easy to state: they fail to relate the even cases 
{ U 2a +4,2b+2 I a > b > 0} to the odd cases { U 2a +3,2b+i | a > b > 0}. 

To remedy this failure, it was sufficient to derive the further^ relation 

n h\ i ft u\ srt , h im /Ca(o + 6- s,s) ( B (a + b- s,s)\ 

C (M)+C(-a,-&) = E(a + 6- S -l)!| (a _ g)!(& _ l)! + (b _ g)!(a _ ^ j (32) 

Ca(M) = aa,b)+a-a,b)-2 1 - a {((a,b) + ((a + b)} (33) 
CbM) = 2 1 " a C(a,6) (34) 



with a > 1 and b > 1. To prove fl32|) , one writes £(a, b) + C( — a ; — b) = J2m n 2/(2m + n) a n 



Then (plf), with A = 2m + n and B = n, yields (B3,p4), after some rearrangements. 



5.2 Double sums in knot theory and field theory 



Adjoining ( p2[ ) to relations in M, it was possible to use REDUCE to derive expressions 
for the 3698 double sums up to level 44, in terms of the 121 irreducible double sums 

{ U 2a+3t2b+1 | a > b > 0, a + b < 20} . (35) 



3 Jon Borwein later told me that (|3^,|33],[34|) were known to, though not used by, the authors of 
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The family of positive knots || that gave rise to this investigation has braidQ words 



{a 2 al a+1 a 2 al b+1 \ a > b > 1} , (36) 

whose enumeration satisfyingly matches that of the irreducibles ( |35|) . 

Note that one omits the knots {a 2 af a+1 a 2 a 1 \ a > 0} from the tally of 3-braids 
in (|36|), since Reidermeister moves transform them to {<T 2a+3 \ a > 0}, which are the 
2-braid torus knots, corresponding [|] to the depth-1 irreducibles {C(2a + 3) | a > 0}, 
whose occurrence has been studied to all [^, [56J orders in quantum field theory. 

Correspondingly, the Euler sums { U 2a +3,i | a > 0} do not occur in countertermsQ though 
they may appear in the finite parts of integrals obtained from Feynman diagrams, and 
hence in the relationships between physical quantities, such as the charge and magnetic 
moment of the electron. 

In fact, two of the most impressive perturbative calculations fl3], Q of physical quan- 
tities in quantum field theory produce the polylogarithm Li 4 (l), with the precise com- 
bination of (ln2) 4 and 7r 2 (ln2) 2 terms given in (§) for = | C(4) — 2a(4). The p- 



parameter |T3[ of electroweak theory entails, at three loops, the level-4 term of 

5 4 = -{8f/ 3 ,i + !C(4)} + 0(e) (37) 

in 4 — 2e spacetime dimensions. The level-4 terms in the three-loop contributions to the 
anomalous magnetic moment of the electron, \{g — 2) e , are obtained from [14| as 
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^3,i + f C(4)}(e/27r) 6 , (3* 



where — e is the electron's charge, in units of (eohc) 1 ^ 2 . Thus the lowest-level double-sum 
irreducible, U^i, is prominent in quantum field theory, though absent from counterterms. 

Having seen the importance of double Euler sums in quantum field theory, and their 
relation to knot theory, one should move on to 3-fold alternating sums, since these too 
occur in field theory, as will be shown in Section 7. Unfortunately, the tools are not yet 
available to derive analytically all the relations implied by (|20"D , with k > 2. One must, 
therefore, make do with analytical derivations of most of them. 



5.3 Proven bounds 

At depths 3 and 4, the rigorous bounds (|28|) were proven by implementing the permu- 
tation (RH) and partial-fraction (|3lD procedures in REDUCE, using its solve command. 



The results are conveniently returned in terms of ARBCOMPLEX |22j variables equal in 
number to the undetermined sums. One may then use the output to form a proven, over- 
complete, demon-free basis. This was achieved for all triple sums up to level 14, and all 
quadruple sums up to level 9, which is no small undertaking, as may be judged from the 
facts that two days of CPUtime were insufficient to solve for the 676 triple sums at level 
15, and 128 megabytes of core memory were insufficient to solve for the 1120 quadruple 
sums at level 10. I recommend these sectors as test grounds for improved algorithms. 

6 For an introduction to knot theory, try |j(J, followed by |3^ |. 

7 For details of renormalization procedures, see fl37|, pending publication of 113] . 
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In the case of triple sums up to level 14, overcomplete bases were constructed by 
adjoining the non-alternating irreducibles {C(5, 3, 3), £(7, 3, 3), ((5, 5, 3)} to the set 

{({2n + 3, —2m — 1, — 2p — 1) \ n + m + p<4, min(n, m, p) > 0} , (39) 

which is also undetermined by the permutation and partial-fraction identities. Attempts 
to reduce this set analytically, by adding identities obtained in the manner of (p2|) , were 
not successful. Therefore the most pressing challenge is to prove the MPPSLQ result that 
E(7, 3) < 2, since only the bound E(7, 3) < 3 has been obtained rigorously, from fl3"9"|). In 
the case of quadruple sums up to level 9, the analytical method fails to find two MPPSLQ 
relations in 5s,4, one of which is given in (P9"|). Hence I recommend the SV,3 and 5*8,4 
sectors as places to start the hunt for more powerful analytical techniques. 



6 Concrete bases 



To aid the elucidation of (PJlq) in Section 8.3, it is instructive to summarize the results of 
MPPSLQ, from Section 4, and REDUCE, from Section 5, by giving concrete irreducibles 
whose values and products span, inter alia, all those spaces of Table 3 with S(l, k) < 36. 
Also included is a choice of the 5 irreducibles in 510,4, which proved attainable with 
MPPSLQ, despite the large size, 5(10,4) = 36. To save space, only argument strings are 
given, with a bar denoting a negative argument, and hence an alternation of sign at the 
corresponding layer of the nest. 

For k = 1, one obviously needs In 2, it 2 , and the odd-zetas, with argument strings 
{(2n+l) | n > 0}. For k = 2, the argument strings {{2n + 1, 2m + 1) | n > m > 0} suffice, 
to level I = 44, and presumably for ever. For the remaining spaces with 5(Z, k) < 36, see 
Table 4, which also includes 5io,4- 

Table 4: Concrete sets of argument strings, yielding demon-free minimal bases. 



55,3 


(3,1,1) 










57,3 


(5,1,1) 


(3,1,3) 








5g,3 


(7,1,1) 


(5,1,3) 


(3,1,5) 






5ll,3 


(9,1,1) 


(7,1,3) 


(5,1,5) 


(5,3,3) 


(3,1,7) 


5l3,3 


(11,1,1) 


(9,1,3) 


(7,1,5) 


(7,3,3) 


(5,1,7) 


56,4 


(3,1,1,1) 










58,4 


(5,1,1,1) 


(3,3,T,T) 








5l0,4 


(7,1,1,1) 


(3,3,3,1) 


(5,3,1,1) 


(5,1,3,1) 


(3,1,5,1) 


57,5 


(3,1,1,1,1) 










58,6 


(3,1,1,1,1,1) 











(5,3,5) (3,1,9) 



It is emphasized that the choices of Table 4 are not the most efficient for integer- 
relation searches; other choices, as in (|1~7D , result in smaller euclidean norms from suc- 
cessful searches. However, one learns more from the patterns of Table 4 than from the 



economic choices of (17). 

First, note that only odd arguments are required: the pattern begun by odd-zetas at 
depth 1 persists. Secondly, note that all arguments may taken as negative, save the first, 
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when the depth, k, is odd. Since the first argument cannot be unity, demons are thereby 
eliminated from the irreducibles. Each is thus computable, to 800 significant figures, in 
half an hour, using fl2l|) . Finally, note that not all permutations of odd integers occur. 

Whether these features persist in larger search spaces, where MPPSLQ may well need 
more than 800 significant figures, can only be conjectured. The argument of Section 8.3, 
however, lends weight to the belief that the key to understanding irreducibility, at level I 
and depth k, lies in an analysis of a restricted set of permutations of partitions of / into 
precisely k positive odd integers. 



7 Euler sums in quantum electrodynamics 



To exemplify the utility and accuracy of the database of exact result for all sums with 
I < 7, it was used to evaluate the Laurent expansion of the three-loop on-shell charge- 
renormalization constant of dimensionally regularized [39| quantum electrodynamics, 
up to terms of level 7, which may be taken as indicative of the transcendentality content 
of four-loop contributions to the anomalous magnetic moment of the electron. 



Like the three-loop corrections |K| to the p-parameter of electroweak theory, on-shell 



charge renormalization involves ||38 
(/i -!)(//-§ 



7T 3 T 3 (3 - fi) 



(p 2 + l)(q 2 + l)(r 2 + l)(p - q) 2 (p - r) 2 
f 1 1 1 



x 



(p — q — r) 2 + 1 (q — r)' 



(40) 



which is a difference of three-loop massive bubble diagrams, regularized in 2/i = 4 — 2e 
euclidean spacetime dimensions. In [^8j it was reduced to a %F 2 hypergeometric series |4(J : 



7_j^ _ r(l-e) r 2 (l + 2e)r(l + 3£) | 
r 2 (l + e)r(l + 4e) J 

TP (1 1 pi p.. 3 I „ 3.i 

3-^2l- L 52 '2 ' 2 ' 2 ' 



24e 4 
+ 



.2 

3e"2 



rr r(l + 2e)r(l + 3e) 



2 6e r 5 (l + e) 



3e 2 (l + 2e) 



(41) 



Expanding the T-functions in the summand of the series, one obtains the e-expansion of 
B4 from boxed sums of the form 



5'odd(a; &!)•••, &fe-l) — X! 



1 



^ (2n - 1)« 



fc-l n-l 

n e 

j=l ITL 



1 



x (2m 4 - 1) 6 > 



(42) 



which, like (|24|), is symmetric in though now it involves reciprocal powers of odd 
integers. To relate these to ([[]), one combines 2 k fc-fold Euler sums with arguments 
of differing sign, to restrict the nested summations to odd integers. Symmetrizing over 
(A; — 1)! permutations of all but the first argument, one exhausts those summations in (H2f) 
with distinct values of {rrii}. Adding combinations of Euler sums with lesser depth, one 
iteratively includes the degenerate summations. Hence the expansion of B 4 to 0(e 3 ) can 
be achieved in terms of the basis (0), by routine (and repeated) reference to the database. 
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The result is 



B 4 = |{-13C(4) + 32a(4)} + |£{-239C(5) + 192a(5) + 2041n2C(4)} 
+ \e 2 {l3 C(6) - 74 C 2 (3) + 160 U 5>1 + 384 a(6) - 204 In 2 2 C(4)} 

+ j^e 3 {-329385 C(7) + 45853 C(3) C(4) + 10875 C(2) C(5) + 7680 C(5, 1, -1) 
-9280 C(3, 3, -1) - 11520 a(4) C(3) + 176031 In 2 C(6) - 2930 In 2 ( 2 (3) 
-48000 ln2 U 5> i + 248832 a(7) + 44064 In 3 2 C(4)} + 0(e 4 ) . (43) 

Note that the polylogarithms (|i"8|) simplify the expansion; had one merely used the con- 
ventional |Il[] polylogarithms Li n (j^J, there would have been 10 additional terms, to this 
order in the expansion. The terms involving In' -4 2 £(4) follow a(l), being generated by 



£(6£ 

Z>4 



J-4 



16 a(/) 
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In2)'- 4 C(4) 



(/-4)! 



(44) 



A strong check of ([I3J) was performed, using the representation 

\ 1 2 



;i + 2e) 



r 1 



r(l + e)r(l-2e) 



1 f 1 dx dy 



o Jo 1 



x 2 y 2 



(l_ x 2 )(1 _ y2) 2 

(8xy) 2 



(45) 



of the hypergeometric series in (f4T|), expanding the integrand to 0(e 5 ), and using the NAG 
routine D01FCF to evaluate 6 double integrals to 8 significant figures. By contrast, 800 
significant figures were rapidly obtained by expressing sums of the form (|42|) in terms of 
those in the database. The agreement of NAG with the more powerful methods developed 
here gives one confidence in the computer-algebraic book-keeping that produced (f|3|). 



There is a further |42|] non-trivial diagram entailed in three-loop charge renormaliza- 
tion: the two-loop fermion-propagator [13, |44| diagram, with a three-fermion intermediate 
state, contributing to on-shell mass renormalization p3l H5j. However, this is eventually 



expressible |38j in terms of T-functions and the same hypergeometric series as is encoun- 
tered in (|4l|) . Thus no further analysis of irreducible Euler sums is entailed, though 
there is a great deal of book-keeping to perform, to obtain the three-loop terms in the 



dimensionally regularized charge-renormalization constant [38 



oo 

l + E^n 

n=l 



e 2 r(i + £ r 

(47r) 2 - £ m 2e / 



(46) 



where e and m are the on-shell charge and mass. 

The one-loop and two-loop contributions are easily found exactly [f44| : 

4 , . 4(1 + 7e - 4e 3 ) 



Ci = - 



Co 



3e J * e(2-e)(l-2e)(l + 2e) 
The first three [3^] terms in the Laurent expansion 

368 t/ „ s 4867 



(47) 



Co 



9e 2 27e I SV ; 3 SV ; 3 SV ; 81 J ^ ' 



(4J 
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are free of contributions from B4, whose expansion starts at level / = 4. Level-3 terms are 



sufficient for the analysis Eg, 47, 48 of a restricted set of four- loop contributions to the 



anomalous magnetic moment of the muon j|9| , and for the study of quark-mass effects at 



the three-loop level of quantum chromodynamics pD |. Thanks to (0), the expansion (58 
may now be continued to level 7, i.e. to the same level as is reached at four loops in the 
electron's anomaly, whose polylogarithms attain order 2L — 1 at L loops. After much 
computer algebra, I was able to obtain: 

C 3 ,, = -3841n 2 2C(2) + ^C(4)-^a(4) + 10241n2C(2)-^C(3) 
_ 1960 104113 

3 SV ; 486 K J 

« -179.724 615 842 918 120 241 823 332 320 650 562 692 071 547 121 , 

C 3 , 2 = 7681n 3 2C(2) - 2841n2C(4) + 384C(2)C(3) + ^C(5) - 1472a(5) 

3 

- 3072 h, 2 2 <(2) + 5!|1 <(4) - ^ <*(4) + 4096 In 2 C(2) 
_3^7 C(3) _Z|4 C(2)+ ™ m 

« -427. 138 027 736 892 466 683 630 594 488 509 635 499 554 227 264 , 
C 3j3 = -1152 In 4 2 C(2) + 852 In 2 2 C(4) - 1280 In 2 £(2) C(3) + C 2 (3) + 6356 C(6) 

- 8832 a(6) - 3680 U 5jl + QUA In 3 2 C(2) — In 2 C(4) + -y- C(2) C(3) 

591473 >/ . 60976 /K . , , >/ , 1472549 >/J , 
+ C 5 a - 12288 n 2 2 C 2 + C 4 

450388 , >/ , 733013 A/ x 25226 ^ x 33051769 , , 

a 4 + 13696 n2 C 2 C(3) C 2 + 51 

27 v ; sv ; 72 sv ; 3 sv ; 17496 v ; 

« -1371.792 496 978 355 362 371 049 120 514 541 715 942 648 466 341 , 

C 34 = ^ln 5 2C(2)- 17041n 3 2C(4) + 38401n 2 2C(2)C(3) + ^^ln2C 2 (3) 
5 54 

334273, 92000, rr 7360 /JW/ N 68689 > , , _ . 

- -^2~ ^ 2 C(6) + -g- In 2 C/ 5)1 + — a(4) C(3) + C(2) C(5) 

334907 , . , , 2459549 . . . . 14720 , 

~ ~W C(3) C(4) + ~^Q~ C(7) " 52 " 2 a{7) " "9" C(5 ' h _1) 
^3360 

+ —=- C(3, 3, -1) - 9216 In 4 2 C(2) + 34067 In 2 2 C(4) - 10240 In 2 C(2) C(3) 

129179 1^944(1 
+ ^C ! (3) + 33767<(6) - 121952 a(6) - ^f^U SA 

+ 24576 2 C(2) _ ™ ln 2 C(4) + Hip C(2 , C (3) 4- ^ C(5) 

_ a(5) _ 41088 ,„ 2 2 ((2) + ^ f (4) - «^ „(4) 

w s 14076461 w . 82735 w . 555842827 , s 

+ 43712ta2«2)--^C<S)-— CW + ^jj- (52) 

« -3273.919 335 883 520 406 469 573 320 145 714 810 021 184 454 681 . 



Terms (H® 

are in agreement with the polylogarithms of order up to 5 that were 



obtained in ||38|| , though they are much simpler in form, thanks to the use of (|I8D 
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Terms are new, and show how complex perturbation expansions may become, 

when they entail polylogarithms of orders 6 and 7, as undoubtedly happens in the four- 
loop contributions to the electron's anomaly. However, I would be rather surprised were 
the four-loop anomaly to contain further transcendentals, not included above, since Euler 
sums are the natural structures to emerge from e-expansions of generalized hypergeo- 
metric series [fi"0"| , |5"T| , [5*2]] whose parameters differ from half-integers |38|, 53, J5~|], or inte- 
gers H, |56j, by multiples of e. In fact, there is good reason to suppose that not all 



of the terms above will occur in the four-loop anomaly, since Laporta and Remiddi WA 
have shown that Li 5 Q) is absent at three loops. 



8 Consequences 



The enumerations (f^[T(J,[16]) have consequences for field theory, knot theory, and number 



theory, which will be discussed in that order. 



8.1 Field theory 

For calculational quantum field theorists, the enumeration fllQD amounts to 'counting the 
enemy', since each new irreducible Euler sum corresponds to the existence of a class 
of polylogarithmic integrals that cannot be related to previously evaluated integrals, by 
computer-algebraic methods. They are few in number, at the levels where calculations 



are currently performed. Their tally, in Table 2, is the integer sequence M0317 of [IS . 

As reviewed in |57| , [58| , |59|1 , there has been tremendous progress in the use of computer 
algebra, most notably by recursive methods within the framework of dimensional regu- 
larization, which automate the computation of single-scale massless [f3("|, and on-shell 

num- 



massive [[HL pq, Eq] propagator diagrams by symbolic manipulation of vast [52 



bers of polynomials in the analytically continued |3J1 dimensionality, 2fi = 4 — 2e, of 
spacetime. The residual difficulty then resides in extracting the Laurent expansions, as 
e —>■ 0, of a small set of irreducible Feynman integrals. Recent success at three-loops JTJ], 



with the electron's anomalous magnetic moment, provides an impressive example of the 
power of this technique. Following discussions at the AI-HEP-92 workshop H , Laporta 
and Remiddi found it possible to achieve a computer-algebraic reduction of all 3-loop 
electron-anomaly diagrams to merely 18 terms \\L4\ , using integration by parts |60], [BTJ in 



4 — 2e dimensions. Drawing on experience of polylogarithmic integration, from previous 
4-dimensional work, they were able to extract the requisite Laurent expansions as e — > 0. 
Along with the inevitable depth-1 irreducibles, In 2 and { £(n) | n < 5}, their result con- 
tains just one irreducible double Euler sum, namely Uz,x, though it is presented in |L4[| in 
terms of the more conventional f|l| polylogarithm Li^lJ, along with the precise (ln2) 4 
and 7r 2 (ln2) 2 terms of (|6|). 

In the case of field-theory counterterms, up to L loops, the irreducible Euler sums 

{C(2a+ 1) | L -2 > a > 0} , (53) 
{ C^o+1,26+1 |a>6>0, L-3>a + 6}, (54) 
{ C(2a + 1, 26 + 1, 2c + 1) | a > 6 > c> 0, a>c, L-3>a + 6 + c}, (55) 
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appear. The first set, studied in 
with the trefoil knot 3i ~ £(3) first appearing in counterterms at L 



53] , |3"4"| , |3J|, corresponds 



to 2-braid torus knots, 
3 loops. The 

second [[16], [HJ set corresponds || to a restricted set of positive 3-braids, with 819 ~ 
first appearing at L = 6 loops, and the third |2], ^[ to positive 4-braids, with the uniquely 
positive 11-crossing non-torus knot o"iO r 2 cr i <T 3 cr 2 cr 3 — C(5, 3,3) first appearing at L = 7 
loops. Hence the 5-loop renormalization of 4 theory was accomplished |J3] in terms of 
only depth- 1 irreducibles, whereas at 6 and 7 loops one first encounters || 2-fold and 
3-fold Euler sums, respectively. 



It is not yet known whether Euler sums exhaust the transcendentals in counterterms 
at L > 7 loops. In addition to IO124 = o 2 a\a 2 a\ ~ Ur,3, there are two further positive 
knots with 10 crossings, namely IO139 = cr 2 cr i (T 2 "i an d IO152 = cr 2 cr i cr 2 cr i) which 7- 
loop analysis || suggests are not associated with Euler sums of depth k < 4. All those, 
and only those, subdivergence-free 7-loop </> 4 diagrams whose link diagrams skein to IO139 
and IO152 appear, on the basis of numerical evidence, to give counterterms that cannot be 
reduced to Euler sums with depth k < 3. Hence it is an open (and fascinating) question 
whether two0 of the 5 new irreducibles in 5*10,4 are associated with these knots. 

Turning to the finite parts of Feynman diagrams, one learns from three-loop analy- 

~~ |). Emphatically no claim is made that Euler 



scs 



13, 14 



that appears, via 
sums exhaust the transcendentality content of perturbative quantum field theory; poly- 
logarithms of non-trivial mass and momentum ratios are everywhere dense. In single-scale 
process, however, where such ratios are unity, or zero, it may occur that the results entail 
only Euler sums, as in the case of electron's anomaly, \{g — 2) e at three loops |14]]- On the 



other hand, with a different configuration of unit and zero masses the maximum value of 



Clausen's integral [4T| is often generated [[42], [53] . 



The key to deciding whether a result is reducible to Euler sums is an analysis of the 
hypergeometric functions [§, ESI |5B[] p+ iF p (l, {a«}; z) that are produced in 4 — 2e 
dimensions. If z = ±1, and {<2j,&i | i < p} differ from integers, or half integers, by 
multiples of s, then reducibility to Euler sums is guaranteed, as in ( ^T] , ^5| ). Such a 
reduction to hypergeometric series has been achieved for the electron anomaly, \{g — 2) c , 
at two loops |38, 



j4j, for the charge-renormalization constant, Z$, at three loops |38f . 
for the Gell-Mann-Low function of quantum electrodynamics to all orders in the large- 
iV limit [0, and for a corresponding limit of the quantum chromodynamics of heavy- 
quark interactions [Q. That is the basis for my strong belief that a similar reduction 
to hypergeometric series, or some generalization |){| |8|, |69j of them, underlies the 
finding of Laporta and Remiddi that the three-loop anomaly involves just one Euler sum 
with depth k > 1, namely the very specific polylogarithmic combination @. It is also 



reasonable to expect that this reducibility to Euler sums persists beyond three loops, 



where only numerical [70] estimates are currently available. 



8.2 Knot theory 

To plain knot theorists, the preoccupations of a knot /field-theorist || §, may appear 
rather restricted. The correspondence claimed by Dirk Kreimer is between positive knots 
and the numbers appearing in field-theory counterterms. His process of discovery began 

8 Recall that L^i and U^,i are absent from counterterms; no positive knot has crossing number 4 or 6. 
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with the observation || that the removal of sub- divergences from Feynman diagrams, 
construed as a skein relation between link diagrams encoding the momentum flow, yields 
a counterterm that is rational if the skeining results in the unknot, as in the case of ladder 
diagrams, where cancellations of ^-functions have been demonstrated perturbatively 0] 
and non-perturbatively [71 . It continued with the observation that the 2-braid torus 
knot (2L — 3, 2) is produced by skeining the crossed-ladder diagram that generates 
C(2L - 3) at L > 3 loops. 

To continue the correspondence, it was clearly necessary to do two things: to go to 
loop-numbers, L, higher than the then current limit of L = 5 |p3"| , which is how I became 
involved ||, and to enumerate positive knots, which we[] have taught REDUCE [Q to 
do, up to 17 crossings, barring degeneracy of HOMFLY |52| polynomials. Here attention 
is largely restricted to knots with up to 13 crossings, as in Table 5. Knots with more 
crossings will figure in 



m 



Table 5: Positive prime knots related to Euler sums, via field-theory counterterms. 
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Excluding the 2-braid torus knots that correspond to depth-1 Euler sums, the tally of 
positive knots with 8 to 13 crossings is apparent in the knot column of Table 5. It forms 
the sequence 1, 0, 3, 1, 7, 8, beginning at 8 crossings with the 3-braid torus knot 
819. By contrast, the corresponding sequence for all prime knots with 8 to 13 crossings is 
given by M0851 of 



18, E3 as 21, 49, 165, 552, 2176, 9988, 



, whose richness reveals the 
specificity of the preoccupations of knot /number /field theorists, to whom, of course, the 



'Table 5 was obtained in collaboration with Dirk Kreimer, who will give further details in |12| 



22 



sparsity of positive knots is a delight: it is broadly commensurate with the slow growth 
in number of irreducible Euler sums, at corresponding levels. 

There is no method, as yet, to assign numbers to knots, other than by brute- force 
evaluation of counterterms from diagrams whose skeinings produce the knots. The obsta- 
cle to high-precision evaluation of such counterterms, beyond 6 loops, is apparent from 
the fact that the knots IO139 and IO152 result, at 7 loops, from diagrams whose evaluation, 
via Gegenbauer-polynomial techniques [f72 |, entails 7- fold summations weighted by the 



squares [73[ of 6-j symbols, in such a manner as to make the evaluation-time of a trunca- 
tion at N increase at least as fast as iV 4 . This is in stark contrast with the linear growth 
in Section 4.1, whose computational challenge hence pales into insignificance, compared 
with that required for [Q. 

Fig. 1: A diagram from cutting a 9- loop bubble that skeins to 13-crossing 4-braids. 



Fortunately, there is a class of Feynman diagrams that may, with effort, be reduced 
to triple Euler sums, by analytical methods. A relatively simple example is provided by 
the 8-loop two-point diagram of Fig. 1, which is finite in 4 dimensions. With unit exter- 
nal momentum, a massless propagator l/p 2 ine for each line, unit vertices, and euclidean 
measure ti~ 2 J d 4 pi oop for each loop, this diagram was evaluated by REDUCE as 

G(3,2,2) = -405 ^3,7,3 + ^AT 2i9>2 + 4680 iV 2 , 7i2 + ^PC(13)- 20535 C(H) 
+ 6480 C(9) C(3) + 12680 C(9) + 480 C(7) C(5) - 19500 C(7) C 2 (3) 

- 7200 C(7) C(3) - 79380 C(7) - 14700 ( 2 (5) C(3) - 1200 ( 2 (5) 

- 68160 C(5) C 2 (3) + 38880 C(5) C(3) - 11520 C 3 (3) + 30240 ( 2 (3) , (56) 



by means of the master formula || [15|] 



G(a, b, c) 



x 



i + j + k)\ 




A(p, m, n) 



i+j+k+r 



(57) 
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where A(p, m, n) results from angular integrations over Chebyshev polynomials and 
is 1, or 0, according as whether g = (p + m + n + l)/2 is, or is not, an integer satisfying 
g > max(p, m, n). 



Each irreducible 3-fold sum appearing in (|5q) belongs to one of the two-parameter 
families: 

AW,2n+i,2m + i = C(2m + 1, 2n + 1, 2m + 1) - C(2m + 1) C(2m + 1, 2n + 1) 

^ /2n + 2A;V , 
+ E( 2A; )Cp(2n + 2Jfe + l,2m-2Jfe + l,2m + l) 



fc=i 

n-l 



2m 4- 2A;\ 

~E( oT Cp(2m + 2A: + l,2n-2fc + l,2m + l), (58) 
fc =o \ 2k J 

N 2 m,2n+i,2m = C(2m, 2n + 1, 2m) + C(2m) {C(2m, 2n + 1) + C(2m + 2n + 1)} 
+ E(^ 2jfe JCp(2n + 2Jfe + l,2m-2Jfe,2m) 

^ /2m + 2k\ , , . 

+ E 2k+1 J(p{2m + 2k + l,2n-2k,2m), (59) 

with product terms 

C P (a, b, c) = C(a) {2 C(6, c) + C(6 + c)} , (60) 

whose systematic inclusion, with the combinatoric factors in (^,^), removes all trace 
of the non-knot number ir 2 from (0), and likewise from every diagram G(a,b,c) with 
a + b + c < 11, according to similar, but much lengthier, evaluations of two-point functions 
obtained by cutting bubble diagrams with up to 13 loops, corresponding to knots with 
up to 23 crossings. 

In general, subdivergence-free bubble diagrams, with up to L loops, yield the knot- 
numbers { iV 0+ 2,26+i,a+2 l a >0, b > 0, L — 4 > a + 6} as the very specific combina- 
tions fl58|j59| ) of Euler sums with depths k < 3 and levels / < 2L — 3 ||. They are knot- 
numbers, in the sense of 0, |], |5| , because counterterms, from subdivergence-free diagrams 
that skein to produce the corresponding knots, contain these numbers, and products of 
other knot-numbers, corresponding to factor knots. The counterterms do not contain 



the non-knot irreducibles In 2 and 7r 2 . Thus the combinations of Euler sums in (^,^) pro- 
vide log-free and 7r-free bases for search spaces in which to evaluate classes of counterterms 
from diagrams that skein to 4-braids. For example, it was possible to evaluate all 7-loop 
4 counterterms from subdivergence-free diagrams that skein to 4-braids in terms of just 
3 knot-numbers: C(H), C(5R 2 (3), and iV 3 , 5 ,3 + 7C(5) C 2 (3) = C(3, 5, 3) - C(3)C(5,3), 
corresponding to the 2-braid torus knot (11,2), the factor knot 5i x 3i x 3i, and the 
uniquely positive 11-crossing hyperbolic^] knot criCT2 cr i cr 3 cr I cr 3 • The factor knot 3i x 819 is 
not produced by skeining the link diagrams that encode the momentum flow. However, 
every diagram that skeins to cr 2 afcr 1 cr 3 crfcr| also produces the other two knots, (11, 2) and 
5i x 3i x 3i. Until a method is devised to predict the rational coefficients with which 
knot-numbers occur, rather than determining them empirically, as at present, the associ- 
ation ^iO"! "! ^ "! "! — iV 3! 5 )3 can be made only modulo ((H) and £(5) C 2 (3) terms. On 
this understanding, iV 3) 5 i3 appears in the 11-crossing entry of Table 5. 

10 In |J], this knot, 10i3g, an d IO152, were wrongly called satellite knots; all three are, in fact, hyperbolic. 
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In Table 5, is the first irreducible triple sum to appear, at 11 crossings. That is 

because the knot-numbers (^,^) are not all independently irreducible. Up to level 13, 
the following odd-zeta reductions are obtained, by extending the methods of to ( [57]) : 



^2,1,2 = 


|C(5) 




^2,3,2 = 


f C(7) 




#3,1.3 = 


-JC(7) 




#2,5,2 = 


f C(9) + |C 3 (3) 




#3,3,3 = 


K(9)-|C 3 (3) 




#4,1,4 = 


WC(9)-|C 3 (3) 




4iV 3i5> 3-5iV 2i7i2 = 


-i§|iC(H)-88C(5)C 2 (3) 




#4,3,4 — 5 #2,7,2 = 


f C(H)-80C(5)C 2 (3) 




2iV 5ili 5-5iV 2i7i2 = 


-i§|iC(H)-56C(5)C 2 (3) 




32 iV 4)5)4 + 140 iV 3j7i3 - 525 iV 2j9i2 = 


24 f 5 C(13) 12880 C(7)C 2 (3) 


- 8400C 2 (5)C(3) 


64 iV 5 ,3, 5 - 100 ^3,7,3 + 175 iV 2 , 9 , 2 = 


-673 C(13) + 4400 C(7) C 2 (3) - 


f 3440C 2 (5)C(3) 


8iV 6jli6 -12iV3 i7i 3 + 49iV 2i9i2 = 


-115C(13) + 976C(7)C 2 (3) + 


752C 2 (5)C(3).(61) 



At level I = 2L — 3 > 3, corresponding to L > 3 loops, there are L — 3 knot-numbers of 
the form { # a+ i,26+i,a+i | a > 0, b > 0, a + b = L - 3}, while (|5|) gives \(L - 3) 2 /12] - 1 
irreducibles Hence the knot-numbers (^,^) fail to exhaust the irreducibles, for 
L > 16 loops. Up to L = 13 loops, their sufficiency has been proven, using REDUCE. 

Two independent irreducibles, chosen to be iV 3j7i3 and A^ 5j3) 5, appear in the 13-crossing 
part of Table 5. They are associated with the first two braid words, on the basis of 
intensive skeining of the link diagrams that encode momentum flows in bubble diagrams. 
A hint of the pen-and-paper labour undertaken by Dirk Kreimer is given by his drawings 
in H, which refer to much simpler Feynman diagrams. Since at least two 13-crossing 
knots of Table 5 emerge from skeining any diagram that yields irreducible triple sums of 
level 13, we cannot yet determine which rational combination of the two irreducibles is 
associated with a given knot. By a combination of skeining and inspection of factorizations 
of Alexander |H| polynomials, we arrive at the association of the first two 13-crossing 
knots of Table 5 with, as yet undetermined, linear combinations of N^j^ and N^^. The 
remaining 6 knot-numbers at level 13 cannot all come from 5^3, since Table 1 reveals that 
12(13, 3) = 7, of which the two non-alternating irreducibles have already been accounted 
for. It is possible that further 13-crossing knot-numbers come from S^s, whose size, 
£(13, 5) = 133, puts it out of the reach of MPPSLQ for the foreseeable future. 

At even levels, the associations of Table 5 are made with combinations 

N a , b = C(-a, 6) - CH>, a) , (62) 

of alternating double Euler sums. The very simple constructs N 5t3 and A^ 73 remove 
{ii 2n | n < 5} from diagrams with up to 7 loops that skein to 819 and IO124. Modulo terms 
corresponding to factor knots, these two knot-numbers agree with the findings of f| , where 
results were written, equivalently, in terms of 29 £(8) — 12 £(5, 3) and 94 £(7, 3) — 793 C(10). 
Previously, the numbers in these combinations appeared gratuitous; now they are seen as 
consequences of ignoring the wider world of alternating sums. 
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Calculations at higher loops reveal that A^2 a +5,3 is a knot-number at level 2a + 8. 
However, a multiple of vr 2 ( a + fe + 6 ) must be subtracted from A^2 a +7,26+5 to obtain a knot- 
number for 3-braids with 12 or more crossings. Table 5 shows that this subtraction has a 
very simple form at 12 crossings. The discovery of the general form of the subtraction is 
frustrated by the fact that counterterms associated with even-crossing knots are generally 
much more difficult to calculate than those associated with odd numbers of crossings, 
as witnessed by the table of results in ||, where 10i 39 and 10i 52 conspired to frustrate 
the precise evaluation of all subdivergence-free contributions to the 7-loop /3-function of 
-theory. As these knots entail computation times scaling as N 4 for truncation at N, 
they leave us with a 7-loop result that is known to 'only' 11 significant figures, after 10 3 
CPUhours. By way of a hard-won, but relatively simple, exact numerical result at 8 loops, 

M(2, 2, 1, 1) = 2 14 • 3 iV 9i3 + 2 4 • 3 ■ 5 3 • 7 C(3) C(9) - 2 5 • 3 3 ( 4 (3) - 2 3 ■ 3 2 • 577 C(5) C(7) (63) 

is offered, as the finding of MPPSLQ for the Feynman diagram of Fig. 2. 

Fig. 2: A diagram from cutting an 8-loop bubble that skeins to 12-crossing 3-braids. 



The evaluation ( [63] ) was accomplished by using 4-dimensional Chebyshev-polynomial 

(64) 



expansions, derived in |36j, for the iteratively-defined coordinate-space constructs 

d 4 z 



P n +i(x,y) 



7T 2 Z 2 



P„(x,z) P (z,y); P (x,y) = l/{x-y) , 



which are then combined by REDUCE to perform a radial integration in 

f d 4 y 

M(a 1 ,a 2 ,a 3 ,a i ) = x 4 J — Y[P ai (x,y) = ^A^rii}) R{{ai}; {^}) 



(65) 



giving 1 15] a 4- fold sum, whose radial term, R, is a huge expression involving inverse 
powers of n« and h = ~ n,, while the angular term, A, vanishes unless h is an integer 
greater than any n.j, in which case A is the smallest of the 8 values rij and h — ni. Taking 
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the REDUCE result, TRANSMP produces MPFUN code that implements the brute- 
force accelerator (f2~5|), with demon-number d = 2, yielding enough significant figures for 
MPPSLQ to discover fl63|) . It would be very hard to obtain such results analytically. As 
so often, in field theory, the whole is much simpler than the parts: thousands of 4-fold 
non-Euler sums produce a one-line, 7r-free result, in terms of the alternating double-sum 
knot-number N 9 ^, which cannot be expressed in terms of non- alternating sums. This led 
me to seek and find such things as (^6|), which in turn led to the enumeration ([|). 

Clearly there is a pressing need for a more developed knot /number/field theory, which 
might tell one which Euler (or other) sums in counterterms to associate with which knots, 
without need of laborious calculations of Feynman diagrams. In particular I would dearly 
like to know whether the two undetermined knot-numbers at level 10 are irreducible 
alternating 4- fold Euler sums, residing in Sio,4. Analytical assistance is ardently sought; 
without it, additional numerical work is likely to add little understanding. 

Eventually, it may prove possible to relate Euler sums to positive knots, directly, in 
a way that is consistent with the field-theory route, yet does not oblige one to follow it. 
That is, I suggest, a substantial task, since it is notoriously difficult to derive non-trivial 
statements that apply to all members of a well-defined class of knot. For example, the 
enumeration, by REDUCE, of positive knots up to 17 crossings, is insecure against the 
possibility that two distinct positive knots might have the same HOMFLY polynomial, 
though no example of positive mutation has come to light. 

It may be that the calculational complexities of field-theory counterterms, and the clas- 
sificational complexities of knot theory, are mutual echoes, with which the now-apparent 
combinatoric simplicity of the enumeration of irreducible Euler sums eventually fails to 
resonate. For the present, however, the following correspondences, from computations to 
13 loops, are rather impressive, to my mind. 



Barring cancellations between diagrams, associated with dynamical symmetries, 
such as local gauge invariance |J or supersymmetry ||, the Euler sum ((2L — 3), 
corresponding [|J to the 2-braid torus knot (2L — 3, 2), first appears in anomalous di- 
mensions at L loops. No other irreducible results from sub divergence-free diagrams 
with less than 6 loops, because no other positive knot has less than 8 crossings. 

At 6 loops, 3-braids start to appear. The first of these is 8i 9 ~ iV 5 3. The irreducibil- 
ity of its knot-number was confirmed in [0, with no knowledge of prior developments 
in field theory flIE , |TT . 



At 7 loops, 4-braids start to appear. The first of these is o r iO r 2 <T i cr 3 cr 2 (T 3 — ^3,5,3- 
The irreducibility of its knot-number was confirmed in f2|, following communication 
of its appearance in field theory ||. The 3-braid 10i 24 ~ N 7j3 is also encountered 
at 7 loops ||, in accord with the tally of |J. 

At 8 loops, there appear: a pair of 13-crossing 4-braids, with knot-numbers iV 3 7 3 
and A^5 i3i 5, in accord with the tally of J2|; and a pair of 12-crossing 3-braids, with 
knot-numbers Ng^ and N? t 5 modulo 7r 12 , in excess of the tally of non-alternating 
double sums in [[]]]. The latter pair led, via fl63|) , to (p6|), which resolved the apparent 
conflict between knot /field theory and number theory and then produced the enu- 
merations as contributions from mathematical physics to pure mathematics. 
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5. 



Results up to 13 loops confirm the association of the knot-numbers (^8|,^) with 
4-braids, up to 23 crossing. 



Results on 14-crossing knots, appearing at 9 loops, will be given in detail in JT2 . 
They confirm the appearance of the expected pair of double-sum irreducibles, A^ 113 
and s modulo 7r 14 . For the first time, a truly irreducible 4-fold Euler sum is 
obtained from a Feynman diagram. The associated knot-number is £(5,3,3,3) + 
C(3, 5, 3, 3) - C(3) C(5, 3, 3) modulo vr 14 . 



8.3 Number theory 

Mathematicians have only recently, it appears, made significant extensions of Euler's 
original study of double sums |J74] . Alternating double sums, familiar in field theory [j3"7 



since [|L6| , were studied in [p]]; triple sums, encountered in [|L7| , were studied in @, |75fl ; 



generic non-alternating sums were studied in |E7J, where the sum rule (E3) was obtained 



The extension of these studies into the entire domain of /c-fold Euler sums, at all levels 
I, with all possible alternations of sign, was undertaken, in this work, in an unashamedly 
experimental^ manner, stemming from an urgent need further to develop the connection 
between knot theory and quantum field theory. 

Mathematics, at its purest, relies on little more than the fertile invention of the human 
mind. Mathematical physics often spawns structures of even greater beauty, thanks to 
what Einstein called the 'incomprehensible' comprehensibility of the natural world. As a 
new variation on this oft-repeated theme, field theory has led, from the observation || 
of the rationality of ladder- diagram counterterms, and the skeining of zeta-rich crossed- 
ladder diagrams @], to a uniquely positive hyperbolic 11-crossing knot in 7-loop || coun- 
terterms, and in this work to the Feynman diagram of Fig. 2, whose evaluation (|63"D then 
resulted in the purely mathematical discovery that non- alternating Euler sums require 
alternating sums for their reduction, as witnessed by (|26|). Setting field theory aside, for a 
brief while, the larger universe of alternating Euler sums proved much easier to enumerate 
than its non-alternating restriction, as witnessed by (|^). 

The validity of the enumeration (|9|) is provisional: a year of work and 10 3 CPUhours, 
of the most exhaustive tests of which I and the engines at my disposal are capable, fail to 
reveal the slightest flaw in it. I invite colleagues with larger numerical appetites to test 
it further, in the lively expectation that it will survive. 

What is needed now is the closest thing to proof for which it is reasonable to hope: 
the establishment by deductive methods of the validity of (Q) as an upper bound on the 
number of irreducible fc-fold Euler sums at level /. The more ambitious aim of proving it 
as an identity is unrealistic, until someone develops the machinery for proving, inter alia, 
the irrationality of £ 2 (137)/C(274), anc i a denumerable infinity of suchlike things. The 
more modest proposal of proving that the number of irreducibles is no more than that 
given by (|9]) seems eminently realistic. To anyone disposed to undertake it, I offer the 
following informal restatement: 



n See |7q| for a suggestion of a working definition of experimental mathematics. 
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1. For concrete values of I and k, such that I + k is even, and I > k > 1, form all the 
partitions of I into precisely k odd integers. 

2. For each partition pi, count the distinguishable permutations of these odd integers 
and denote the answer by Pj. 

3. Let Ai be the number of products of lower-level irreducibles, associated with parti- 
tion pi by the generator fll~2]) operating on already established irreducibles. 

4. Summing Ei = Pi — A^ over the odd partitions one arrives at ([]). 

By way of example, consider Sio,4, which admits of the odd partitions pi = 7+1 + 1 + 1, 
p 2 = 3 + 3 + 3 + 1, and p 3 = 5 + 3 + 1 + 1, with P 1 = P 2 = 4!/3! = 4 and P 3 = 4!/2! = 12. 
To exemplify the structure, only the argument strings of Section 6 are notated, with (1) 
standing for — In 2. Thus the 3 products 

(7) (I) (I) (I), (7, 1) (I) (I), (7,1,I)(I), (66) 
leave E\ = 4 — 3 = 1 as the number of irreducibles associated with p%. Similarly, 

(3)(3)(3)(T), (3,T)(3)(3), (3, 1, 3) (3), (67) 

leave E2 = 4 — 3 = 1 irreducibles. Finally, the 9 products 

(5)(3)_(T)(T), (5,3)(T)(T), (5, T)(3)(I), (5 I)(5)(I), (3,T)(5,T), 

(5, 1, 1)(3), (3, 1, 1)(5), (5, 1, 3)(1), (3, 1, 5)(1), (68) 

leave E 3 = 12 - 9 = 3 irreducibles. The tally £(10,4) = 1 + 1 + 3 = 5 agrees with (g). 
Note that in (|66| , |67| , |68D the specific choices of Section 6 were made, for irreducibles in 
spaces of lesser depth and level. However, that was not necessary: all one needs is the 
number of irreducibles, associated with the sub-partition in the smaller space, which is 
itself determined by the Aufbau. By the same token, all one needs to carry forward 
from Sio^, for later stages of the Aufbau, are the numbers, Ei, of irreducibles associated 
with each of the 3 partitions of 10 into 4 odd integers; it is not necessary to specify a 
concrete basis by choosing signs or specific orderings of arguments. Note also that the 

/ l + k-2 \ 

total number of permutations is merely an element of Pascal's triangle: J2i Pi — ( uli ) > 

in general, with 4 + 4 + 12 = (fj in 510,4. It is, therefore, an elementary exercise to build 
up Euler's triangle of Table 1, by repeated application of Ei = Pi — Ai, which is, so to 
speak, the microscopic version of the macroscopic results (P|,|10|)[16D . To paraphrase: 



The number of irreducibles is given, by Euler's triangle, as the deficit between 
the number that Pascal's triangle generates, by permutation of integers in odd 
partitions, and the number of products already given by previous deficits. 

In its final form, a proof of @, as a rigorous upper bound, may appear almost as simple 
as the "Euler = Pascal - Already" paraphrase of what is to be proved. It may, however, 
require formidable combinatoric intuition to generate and, more importantly, to organize 
sufficient permutation and partial-fraction relations, between /c-fold Euler sums, with all 
possible alternations of sign, as effectively as the authors of organized the relations 
between non-alternating sums at depths k = 2 and k = 3. Mindful of what was entailed 
by these lesser tasks, I did not attempt the greater. 
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8.4 Conclusion 



I marvel that the quantum electrodynamics [77] of Dyson, Feynman, Schwinger, and 
Tomonaga |78j leads, after 50 years of dedicated calculation and the attainment of equally 



impressive experimental accuracy [f79[ , to agreement [£L4[ between theory and experiment 
on the eleventh significant figure in the magnetic moment of the electron. That skillful 

would smile, one feels, on seeing nothing 
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calculator of double sums, Leonhard Euler 
more complicated than J2n>m{~ l) n+m /n 3 m in the sixth-order perturbation expansion. 
Gauss, too, might be amused to see numbers, from his hypergeometric series f3~8| , |40| , 
and knots, whose codification |K| he began, walk hand in hand, down the perturbation 
expansions of quantum field theory, to all || orders in the coupling constant. 

Such reassuring order, in the mathematical description of nature, at astonishing levels 



of accuracy, reinforced by experience M, I I 



3], |54| with hypergeometric 



series generated by dimensional regularization, and now by the remarkable matches be- 
tween Euler sums |jl], § and Dirk Kreimer's knot /number/field theory || |], [5], ||, || , leads 
me to believe that the k < 3 sums in and hence in (0), more than suffice for the 
electron's anomaly at eighth order. Certainly, they suffice for the reduction of all Euler 
sums with levels / < 7. At any level, one has only to consult (|9]) for the tally of irreducible 
/c-fold Euler sums, be they construed as calculational obstacles or mathematical friends. 
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